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PACS 89 . 75 . Da - Systems obeying scaling laws 
PACS 05 . 45 . Df - Fractals 

PACS 02. 10. Ox - Combinatorics; graph theory 

PACS 89 . 75 . Fb - Structures and organization in complex systems 

Abstract. - Many real networks share three generic properties: they are scale-free, display a 
small-world effect, and show a power-law strength-degree correlation. In this paper, we propose a 
type of deterministically growing networks called Sierpinski networks, which are induced by the 
famous Sierpinski fractals and constructed in a simple iterative way. We derive analytical ex- 
pressions for degree distribution, strength distribution, clustering coefficient, and strength-degree 
correlation, which agree well with the characterizations of various real-life networks. Moreover, 
we show that the introduced Sierpinski networks are maximal planar graphs. 



Introduction. — In the last few years, research of 
complex networks have become a focus of attention from 
the scientific community [1-5]. One of the main reasons 
behind the popularity of complex networks is their flexi- 
bility and generality for representing real systems in na- 
ture and society. Researchers have done a lot of empirical 
Studies, uncovering that various real-life networks shar- 
ing some generic properties: power-law degree distribu- 
tion [6], small- world effect including small average path 
length (APL) and high clustering coefficient [7] . Recently, 
many authors have described some real-world systems in 
terms of weighted networks, where an interesting empirical 
phenomenon has been observed that there exists a power- 
law scaling relation between the strength s and degree k 
of nodes, i.e. s ^ with f3 > 1 [8-11]. 

With the intention of studying the above properties of 
real-world systems, a wide variety of models have been 
proposed [1-4]. Watts and Strogatz, in their pioneering 
paper, introduced the famous small-world network model 
(WS model) [7] , which exhibits small APL and high clus- 
tering coefficient. Another well-known model is Barabasi 
and Albert's scale-free network model (BA model) [6], 
which has a degree distribution of power-law form. How- 
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ever, in these two elegant models, scale-free feature and 
high clustering are exclusive. Driven by the two semi- 
nal papers [6,7], a considerable number of other models 
have been developed that may represent processes more 
realistically taking place in real- world networks [12-18]. 
Very recently, Barrat, Barthelemy, and Vespignani have 
introduced a model (BBV) for the growth of weighted net- 
works [19,20], which is the first weighted network model 
that yields a scale-free behavior for strength and degree 
distributions. Enlightened by BBV's remarkable work, 
various weighted network models have been proposed to 
explain the properties found in real systems [21-27] . These 
models may give some insight into the realities. Par- 
ticulary, some of them presente all the above-mentioned 
three characteristics such as power-law degree distribu- 
tion, small-world effect, and power-law strength-degree 
correlation [21-23]. Although great progresses have been 
made in the research of network topology, modeling com- 
plex networks with general structural properties is still of 
current interest. 

On the other hand, fractals are an important tool for the 
investigation of physical phenomena [28] . They were used 
to describe physical characteristics of things in nature and 
life systems such as clouds, trees, mountains, rivers, coast- 
lines, waves on a lake, bronchi, and the human circulatory 
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system, to mention but a few. A vast literature on the 
theory and apphcation of fractals has appeared. Among 
many deterministic and statistical fractals, the Sierpinski 
gasket [29,30] is one of the earliest deterministic fractals; 
it has provided a rich source for examples of fractal be- 
havior [28,31,32]. Our initial physical motivation for this 
work lies in the use of the Sierpinski gasket as models for 
complex networks. 

In this letter, based on the well-known Sierpinski fam- 
ily fractals, we introduce a class of deterministic networks, 
named Sierpinski networks. We propose a minimal itera- 
tive algorithm for constructing the networks and studying 
their structural properties. The networks are maximal 
planar graphs, show scale- free distributions of degree and 
strength, exhibit small- world effect, and display power-law 
strength-degree correlation, which may provide valuable 
insights into the real-life systems. 

The network derived from Sierpinski gasket. — 

We first introduce a family of fractals, called Sierpinski 
fractals, by generalizing the construction of the Sierpinski 
gasket. The classic Sierpinski gasket, shown in Fig 1(a), 
is constructed as follows [29]. We start with an equilat- 
eral triangle, and we denote this initial configuration by 
generation t = 0. Then in the first generation t = 1, 
the three sides of the equilateral triangle are bisected and 
the central triangle removed. This forms three copies of 
the original triangle, and the procedure is repeated indef- 
initely for all the new copies. In the limit of infinite t 
generations, we obtain the well-known Sierpinski gasket 
denoted by SG2{t). Another fractal based on the equilat- 
eral triangle can be obtained if we perform a trisection of 
the sides and remove the three down pointing triangles, to 
form six copies of the original triangle. Continue this pro- 
cedure in each copy recursively to form a gasket, denoted 
SG3{t), shown in Fig. 1(b). Indeed, this can be general- 
ized to SG^^{t), for any positive integer w with w ^ 1, by 
dividing the sides in lu, joining these points and removing 
all the downward pointing triangles [33]. Thus, we obtain 
a family of fractals (Sierpinski fractals), whose Hausdorff 
dimension is dj — log[|a;(a; -|- l)]/log(a;), which tends to 
2 as w ^ oo [34]. 

One can use these fractals to construct networks. The 
translation from the fractals to network generation is quite 
straightforward. Let the nodes (vertices) of the networks 
correspond to the removed triangles and make two nodes 
connected if the boundaries of the corresponding triangles 
contact one another. For uniformity, the three sides of 
the initial equilateral triangle at step also correspond to 
three different nodes. Figure 2 shows a network based on 
SG3{2). All generated networks have similar properties. 
In what follows we will focus on the particular network 
corresponding to the case of = 3 (sec Figure 2), which 
is also called Sierpinski network, and other cases can be 
treated analogously. 

Iterative algorithm for the maximal planar net- 
work. — In the construction process of SG3{t), for each 
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Fig. 1: The first two stages in the construction of two fractals 
from the family of Sierpinski gaskets. 

equilateral triangle at arbitrary generation, once we per- 
form a trisec;tion of its sides and remove the three resultant 
down pointing triangles, six copies of it are formed. When 
building the network, it is equivalent that for every group 
of three new added nodes, six new small equilateral trian- 
gles are generated, each of which may create three nodes in 
the next generation. According to this, we can introduce 
a general algorithm to create the corresponding network, 
denoted by F{t) after t generation evolutions. 

The iterative algorithm for the network is as follows: 
For t = 0, F{0) consists of three nodes forming a triangle. 
Then, we add three nodes into the original triangle. These 
three new nodes are linked to each other shaping a new 
triangle, and both ends of each edge of the new triangle 
are connected to a node of the original triangle. Thus we 
get F{1), see Figure 3. For t > 1, F(t) is obtained from 
F{t—1). For each of the existing triangles of F (t — 1 ) that 
is not composed of three simultaneously emerging nodes 
and has never generated a node before, we call it an active 
triangle. We replace each of the existing active triangles 
of F{t — 1) by the connected cluster on the right hand of 
Figure 3 to obtain F{t). The growing process is repeated 
until the network reaches a desired order (node number of 
network). Figure 2 shows the network growing process for 
the first two steps. 

Next we compute the order and size (number of all 
edges) of the network F{t). Let Ly{t), Le{t) and L^{f?) be 
the number of vertices, edges and active triangles created 
at step t, respectively. By construction (see also Figure 3), 
each active triangle in F{t — 1) will be replaced by six ac- 
tive triangles in F{t). Thus, it is not difficult to find the 
following relation: L^it) = 6L/^{t — 1). Since L^{0) = 1, 
we have L/\{t) ~ 6*. 

Note that each active triangle in i^(i — 1) will lead to an 
addition of three new nodes and nine new edges at step t, 
then one can easily obtain the following relations: Ly{t) = 
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Fig. 2: Illustration of a deterministic network in relation to 
SG3(2). 



3LA{t-l) = 3 -6*-!, and Le{t) =9LA{t-l) = 9 • 6*-^ 
for arbitrary t > 0. From these results, we can compute 
the order and size of the network. The total number of 
vertices Nt and edges Et present at step t is 



3 • 6* 12 



ti=0 



and 



Et = ^L,{U) = 



9 • 6* -f 6 



(1) 



(2) 



2Et 
Nt 



respectively. So for large t, the average degree fct 
is approximately 6, which shows the network is sparse as 
most real systems. 

From Eqs. (1) and (2), we have Et = 3Nt - 6. In addi- 
tion, by the very construction of the network, it is obvious 
that arbitrary two edges in the network never cross each 
other. Thus our network is a maximal planar network (or 
graph) [35], which is similar to some previously studied 
networks [36-39]. 

Relevant characteristics of the network. Now 

we study the statistical properties of the network, in terms 
of degree distribution, clustering coefhcient, average path 
length, and strength distribution. 

Degree distribution. When a new node i is added to 
the graph at step ti {ti > 1), it has a degree of 4. Let 
L^iijt) be the number of active triangles at step t that 
will create new nodes connected to the node i at step t + 
1. Then at step ti, L/\{i,ti) = 3. From the iterative 
generation process of the network, one can see that at 
any step each two new neighbors of i generate three new 
active triangles with involving i, and one of its existing 
active triangle is deactivated simultaneously. We define 
ki{t) as the degree of node i at time t, then the relation 



Fig. 3: Iterative construction method for the network. 



between ki{t) and L/^li^t) satisfies: 



iA(i,<) = h{t) - 1. 



(3) 



Now we compute LA{i,t). By construction, LA(i,t) — 
3 L^ii, t—1). Considering the initial condition L^ii, ti) — 
3, we can derive L^iiji) = 3*^**"'"^. Then at time t, the 
degree of vertex i becomes 



fc^(t) =3 



t-ti+i 



1. 



(4) 



It should be mentioned that the initial three vertices 
created at step have a little different evolution process 
from other ones. We can easily obtain: LA{0,t) — 3* and 
ki{t) = 3' + 1. So at step t, initial three vertices have the 
same degrees as those ones born at step 1. 

Equation (4) shows that the degree spectrum of the net- 
work is discrete. It follows that the cumulative degree 
distribution [3] is given by 



Pr, 



L,{t) 3-6* 



12 



T<ti 



Nt 



3 • 6* + 12 



(5) 



Substituting for ti in this expression using ti = t + 1 — 

In(fc-l) . 

gives 

18.6*.(fc~l)-('"^/'"^) + 12 
PcumW - 3 . 6* + 12 • 

When t is large enough, one can obtain 

Pcu„.(fc) = 6.(fc-l)-[i+(i-2/'"3)l. (7) 

So the degree distribution follows a power law form with 
the exponent 7^ = 2 -f . 

Clustering coefficient. The clustering coefficient [7] Ci 
of node i is defined as the ratio between the number of 
edges Ci that actually exist among the ki neighbors of node 
i and its maximum possible value, ki{ki — l)/2, i.e., Ci = 
2ei/[ki{ki — 1)]. The clustering coefficient of the whole 
network is the average of CjS over all nodes in the network. 
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Fig. 4: Semilogarithmic plot of average clustering coefEcient 
Ct versus network order Nt. 



Fig. 5: Average path length dt versus network order Nt on a 
semilogarithmic scale. The solid line is a guide to the eye. 



For our network, the analytical expression of clustering 
coefficient C{k) for a single node with degree k can be 
derived exactly. When a node enters the system, both fc, 
and Ci are 4. In the following iterations, each of its active 
triangles increases both fcj and Cj by 2 and 3, respectively. 
Thus, Ci equals to 4 + | (fci — 4) for all nodes at all steps. 
So one can sec that there exists a one-to-one correspon- 
dence between the degree of a node and its clustering. For 
a node of degree k, we have 



C{k) = 



2e 



2[4+|(/c- 



k{k-l) 



k{k-l) 



k 



1 



k-1 



(8) 



In the limit of largo k, C{k) is inversely proportional to 
degree k. The same scaling of C(fc) ^ has also been 
observed in several real- life networks [40]. 

Using Eq. (8), we can obtain the clustering Cf of the 
networks at step t: 



r=0 



4 



Dr - 1 



(9) 



where the sum runs over all the nodes and Dr is the degree 
of the nodes created at step r, which is given by Eq. (4). 
In the infinite network order limit {Nt oo), Eq. (9) 
converges to a nonzero value C = 0.598, as shown in Fig. 4. 
Therefore, the average clustering coefficient of the network 
is very high. 

Average path length. Shortest paths play an impor- 
tant role both in the transport and communication within 

a network and in the characterization of the internal struc- 
ture of the network. We represent all the shortest path 

lengths of F{t) as a matrix in which the entry dij is the fj-gni which we have 
geodesic path from node i to node j, whore geodesic path 
is one of the paths connecting two nodes with minimum 
length. The maximum value of dtj is called the diameter of 
the network. A measure of the typical separation between 



two nodes in F{t) is given by the average path length dt, 

also known as characteristic path length, defined as the 
mean of geodesic lengths over all couples of nodes. 

In Fig. 5, we report the dependence relation of APL dt 
on network size N^. From Fig. 5, one can see that the 
average path length grows logarithmically with increasing 
size of the network. This logarithmic scaling of dt with 
network size Nt, together with the large clustering coeffi- 
cient obtained in the preceding subsection, shows that the 
considered graph has a small- world effect. 

Strength distribution and strength- degree correlation. 
Strength usually represents resources or substances allo- 
cated to each node, such as wealth of individuals in finan- 
cial contact networks [41], the number of passengers in 
airports in world-wide airport networks [42] , the through- 
put of power stations in electric power grids [43], and so 
on. In our model, the strength of a node is defined as the 
area of the removed triangle it corresponds to [44]. For 
uniformity, let the initial three nodes born at step have 
the same strength as those created at step 1. 

We assume that the area of the initial equilateral trian- 
gle of the Sierpinski gasket is A. By the very construction 
of the network, all simultaneously emerging nodes have 
the same strength, because their corresponding triangles 
have identical area. It is easy to find that each removed 
triangle covers the portion = i of one removed triangle 
in the preceding generation. After t iterations, all nodes 
which are generated at a certain step ij = 1,2, - ■ ■ ,t have 
the strength Si{t): 



i{t) = M** 



(10) 



Si{t) = /X*'-* • A, (11) 
where A is the area st{t) of a triangle removed at step t. 
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Eqs. (4) and (11) yield a power-law correlation between 
strength and degree of a node: 



If, (12) 



which implies when k is large enough, s ^ k"^. This non- 
trivial power-law scaling between strength of a node and 
its degree has been empirically observed in a variety of real 
networks, such as the airport networks, the shareholder 
networks, and the Internet. 

Analogously to computation of degree distribution, we 
one can find that the strength distribution is also scale- free 
with exponent 7s as: 



7s 



In 2 
21n3" 



(13) 



As known to us all, for weighted networks with non- 
linear strength-degree correlation s ~ k^, if their dis- 
tributions of degree and strength behave as power laws, 
P{k) ~ k~'"' and P{s) ~ s"'''", then there is a general 
relation between 7^ and 7s as [19,20]: 



~ /? ^ /? • 



(14) 



We have shown that in our model 71, = 2 + ^ and B = 2. 

In 3 

According to Eq. (14), the exponent of strength distribu- 
tion is 7s = I + ^f^, giving the same 7s value as that 
obtained in the direct calculation of the strength distribu- 
tion, see equation (13). 

Conclusion. — Deterministic model makes it easier 
to gain a visual understanding of how do different nodes 
relate to each other forming complex networks [45 50]. 
On the basis of Sierpinski fractals, we have proposed and 
studied a kind of deterministic networks. According to 
the network construction processes we have presented an 
algorithm to generate the networks, based on which we 
have obtained the analytical results for degree distribu- 
tion, clustering coefficient, strength distribution, as well as 
strength-degree correlation. We have shown that the net- 
works have three important properties: power-law distri- 
butions, small-world effect, and power-law strength-degree 
correlation, which are in good accordance with a variety of 
real-life networks. In addition, the networks are maximal 
planar graphs, which may be helpful for designing printed 
circuits. 

Although we have studied only a particular network, in 
a similar way, one can easily investigate other Sierpinski 
networks with various values of 7s and 7^ , and their gen- 
eral properties such as small-world effect and power-law 
strength-degree relation are similar. Moreover, using the 
idea presented, one can also establish random networks, 
which display similar features as their deterministic coun- 
terparts studied here. As the classic Sierpinski gaskets are 
important for the understanding of geometrical fractals in 
real systems, we believe that our research could be useful 
in the understanding and modeling of real-world networks. 
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